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Observations of the coronae of the Sun and of solar-like stars provide complementary 
information to advance our understanding of stellar magnetic activity, and of the 
processes leading to the heating of their outer atmospheres. While solar observations 
allow us to study the corona at high spatial and temporal resolution, the study of 
stellar coronae allows us to probe stellar activity over a wide range of ages and 
stellar parameters. Stellar studies therefore provide us with additional tools for 
understanding coronal heating processes, as well as the long-term evolution of solar 
X-ray activity. We discuss how recent studies of stellar magnetic fields and coronae 
contribute to our understanding of the phenomenon of activity and coronal heating 
in late-type stars. 
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1. Introduction 


The detailed physical processes that heat the outer atmosphere of the Sun and 
of solar-like stars to millions of degrees are still poorly understood, and remain a 
major open issue in astrophysics. Recently, significant progress has been made, both 
observationally and theoretically, in understanding coronal heating mechanisms (see 
also other papers of this issue on “New approaches in coronal heating” for reviews 
of different aspects of recent coronal heating research). 

Stellar magnetic activity is ubiquitous in solar-like stars, i.e., relatively un¬ 
evolved stars with convective envelopes below their photospheres, and it includes 
a variety of phenomena such as starspots, activity cycles, hot outer atmospheres 
(comprised of chromosphere, transition region and corona), and magnetic breaking 
due to stellar winds (see e.g., a review by Schrijver fe Zwaan| 2008). Understanding 
magnetic activity is of fundamental importance for stellar astrophysics since it plays 
a crucial role in star and planet formation and evolution, controls stellar angular 
momentum loss, and drives orbital decay and evolution in close binary systems. 

Solar observations provide us with a close up view of magnetic activity, allowing 
us to investigate the characteristics of the magnetic field and magnetized plasma at 
high spatial and temporal resolution. In Fig. [T] we show an example of recent solar 
observations at UV/EUV/X-ray wavelengths, able to resolve coronal structures 
down to ~ 250 km. Stellar observations, on the other hand, allow us to explore how 
the stellar activity depends on stellar parameters (e.g., stellar mass, age, chemical 
composition), and to test the validity and limitations of the solar-stellar analogy, 
and put tighter constraints on models of activity and coronal heating. 
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In this paper we focus on the stellar perspective on magnetic activity, and in 
particular we discuss the insights that can be gained on coronal activity and coronal 
heating from stellar studies. By necessity, we will confine the discussion to a few 
selected topics that we believe have provided more significant advances in recent 
years, and in particular will review results based on observational studies. We re¬ 
fer the reader to earlier reviews for in-depth discussion of other aspects (including 
theoretical modeling) that are not included or only very briefly touched upon here 
(e.g., dynamo models, stellar chromospheric emission and cycles, starspots, activity 
in pre-main sequence stars, star-planet interactions and effects on stellar activity; 


e.g., Durney et al. 1993] Charbonneau[2010| Linsky 1980 Berdyugina 2005; Strass- 

meier 

2009; Favata & Micela|2003| Preibiscli et al. 2005 Giidel & Naze 2009; Testa 

2010& 

Poppenhaeger 20141. 



Figure 1. Example of recent imaging observations of the solar corona taken with the 
Atmospheric Imaging Assembly (AIA-, [Lemen et al. 2012) onboard the Solar Dynamics 


Observatory (left), and with the Interface Region Imaging Spectrograph (De Pontieu et al. 
2014] right). AIA provides full-disk images in narrow EUV passbands, at 12 sec cadence 


and ~ 0.6 arcsec/pix, and IRIS provides simultaneous spectra and images of the photo¬ 
sphere, chromosphere, transition region, and corona, with temporal resolution of seconds, 
~ 0.166 arcsec/pix, and ~ 1 km s -1 velocity resolution, over a field-of-view of up to 175 
arcsec x 175 arcsec. Left.: three-color full disk image combining AIA observations in three 
narrow EUV passbands centered around 94A (red), 193A (green), and 17lA (blue), re¬ 
spectively. Right: slit-jaw images in the IRIS FUV passband centered around 1400A, of 
active region 12192, which is the largest active region observed on the Sun since 1990. 
Both observations are taken on 2014-10-25 around 18UT. 


In Section [2] we briefly review the basic properties of magnetic activity in solar¬ 
like stars, focusing in particular on X-ray observations and coronal heating pro¬ 
cesses. We then discuss a few selected aspects of magnetic activity in some detail, 
including its correlation with stellar parameters and evolution with stellar age. In 
Section [3j we highlight the recent results on activity cycles and magnetic grand 
minima, and in Section [4] we discuss the chemical abundance anomalies in stellar 
coronae and their possible diagnostic potential for coronal heating. Finally in Sec- 
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tion [5j we discuss the results we have highlighted and conclude with a series of key 
open questions. 


2. Stellar magnetic activity, its correlation with stellar 
parameters, and its evolution 

The evidence of a hot outer atmosphere on the Sun was initially interpreted as 
a direct consequence of convection, via dissipation of acoustic modes generated 


by turbulent convection (e.g., Schwarzschild 1948 Kuperus 1969). This scenario 
came into question following the first X-ray observations of the Sun and other 
stars. A rapid improvement of instrumental capabilities since the first X-ray solar 


observations (Friedman et al. 1951) yielded X-ray images of the Sun’s high-energy 
emission. 


and Skylab (e.g., Vaiana et al. 


Blake et al. 

1963 

Vaiana et al. 

1968 


from rocket experiments (e.g, 

1973) These observations 


revealed the highly structured and inhomogeneous nature of the hot corona, and 
their relation to magnetic features (e.g., van Speybroeck et al. |1970). Systematic 


X-ray observations of late-type stars with space observatories (e.g., Vaiana et al. 
1981), revealing the wide range in X-ray luminosity of stars with similar spectral 


types, further demonstrated the inadequacy of acoustic heating in explaining the 
emission from stellar coronae. Instead, magnetic heating seemed to be required (e.g. 
Rosnerl|1980| see also reviews by e.g.JHaisch & Schmitt 1996; Ulmschneider et al. \ 
19911). 


While progress in solar and stellar coronal observations has helped to develop 
and constrain models for coronal heating, the details of the physical processes con¬ 
verting magnetic energy into atmospheric heating are not yet well established. The 
most favored candidate heating mechanisms include dissipation of magnetohydro¬ 


dynamic (Alfven) waves (e.g., van Ballegooijen et al. , 2011 2014), and dissipation of 
magnetic stresses as in the “nanoflare” model, where random photospheric motions 


Parker 

1988 

Priest et al. 

2002) 


2002). Most viable models predict energy release char¬ 


acterized by small spatial and temporal scales, typically below current resolution 
capabilities for the Sun, and constraining the models with observations presents 
significant challenges (see reviews by e.g., Klimchuk 2006; Reale |2014 ). 

The main characteristics of X-ray emission in late-type stars can be summarized 


as follows (see also e.g., Giidel 2004 Testa 2010& for reviews): (1) all unevolved 
late-type stars (F-M type) are X-ray sources (e.g. |Vaiana et al.\ |1981| |Schmitt 
fe Liefkel [2004]), as are evolved late-type stars earlier than spectral type mid-K 


(Haisch et al.\ 1992); (2) they show a very wide range in X-ray emission level of 


several orders of magnitude, even among the same spectral types; (3) the y ap p ear 
to have a minimum X-ray surface flux, Fx, value (~ 10 4 erg enr 2 s _1 ; Schmitt 
& Liefke e.g., 2004); (4) their X-ray activity, as parametrized by the ratio Lx/L boi 


of the X-ray luminosity, Lx, to the total (bolometric) luminosity, L f,oi, appears to 
saturate at a value of log(Lx/Lf, oi ) ~ —3 (Vilhu 1984[ Vilhu & Walter 1987[ see 
sec.2(a] below for a discussion of the implications of some of these properties on 
dynamo action at work in late-type stars). 

As mentioned above, X-ray emission is only one of the aspects of stellar magnetic 
activity, and it is therefore interesting to explore correlations with other activity 
proxies, such as chromospheric activity, and magnetic flux density. Analogous to 


Article submitted to Royal Society 












































































































4 


P. Testa, S.H. Saar, J.J. Drake 


coronal X-ray activity, chromospheric activity also exhibits a large range for stars 
with similar surface gravity and effective temperature (saturation effects for chro¬ 
mospheric emission are discussed in sec.2 ja] below). 

The X-ray flux shows a good correlation with chromospheric activity indicators 


based on Call HK core emission, and with the magnetic flux density (Schrijver 


et al. 1992). These flux-flux relationships hold for stars (early-F to rnid-M type; 


main sequence, subgiants; and giants down to rnid-K), for the Sun as a whole, and 
for specific solar features extending over 2-4 orders of magnitude, depending on the 
diagnostics used (see Schrijver & Zwaan 2008 and references therein). |Schrijver1F 


Zwaan (2008) parametrize the relation between magnetic flux density, and the 
radiative diagnostics of activity, F), as power-laws: F, = ai \4>\ bi ■ They find that 
the exponent bi monotonically increases with the temperature of formation of the 
radiative diagnostics, ranging from ~ 0.5 for chromospheric emission in Can and 
Mg II, to ~ 0.75 for the transition region emission in e.g., SilV and Civ, to ~ 1 for 
the coronal X-ray emission (see also Ayres et al. 1995). 

The correlation of X-ray emission and magnetic flux has been further explored 
by Pevtsov et al. (2003) who find a roughly linear relationship between Lx and the 
total unsigned magnetic flux (Lxoc (J> 113±0 05 ) ) for spatially resolved solar features 
(e.g., active regions, X-ray bright points, quiet Sun regions) and active stars (dwarfs 
and pre-main sequence), over ~ 12 orders of magnitude (see also Fisher et al. 1998). 
These results extend beyond stars for which the coronal filling factor has already 
reached close to unity (i.e., stars that are completely covered in active regions; 
2§, ng that the increase in X-ray emission with magnetic flux is 


see sec 


likely associated with an increase in magnetic field strength and not simply with 
the integrated surface area of emitting regions. Pevtsov et al. (2003) interpret their 
result as being a consequence of a universal relationship between magnetic flux 
and power dissipated through coronal heating, and suggest this might imply the 
presence of one dominant coronal heating mechanism for all solar coronal features 
and stellar coronae (see section [5] for a brief discussion). 

Chromospheric observations of late-type stars also reveal evidence of a “basal” 
emission component that appears unrelated to the stellar magnetic activity, cov¬ 


ering the entire surface of magnetically inactive stars (e.g., Schrijver et al.. 1989 


Schrijver 19951. This basal activity has been interpreted as due to acoustic heating 


(see Schrijver 


1995), however, Bercik et al. (2005) indicate that a turbulent dy¬ 


namo, from non-rotating plasma, represents a viable alternative to acoustic heating 
models, and can reproduce the observed basal coronal (X-ray) and chromospheric 
(Mgn) flux. Cuntz et al. (1999) have demonstrated that MHD wave heating models 
using direct stellar magnetic measurements can successfully heat the chromospheres 
over a wide range of activity. 

The close similarity of the phenomenology of solar and stellar activity consti¬ 
tutes the basis of the “solar-stellar connection” (e.g., Peres et al. 2000 20041. The 
challenge of the solar-stellar connection is to investigate the limitations of the solar 
analogy. In fact, several significant differences in solar and stellar activity phenom¬ 
ena are also observed, especially going towards the extreme of activity. Active stars, 


for example, tend to have large polar spots (e.g., Berdyugina et al. 1998 2000 Hus¬ 


sain et al. 2007), at odds with less active stars like the Sun. Also, the high-level of 


X-ray emission of active coronae, several orders of magnitude larger than that of 
the Sun, cannot be simply reproduced by scaling up the solar corona, i.e., covering 
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the entire stellar surface with solar active regions (e.g., Drake et al. 2000). We 


now know that active stellar coronae differ in at least some aspects from the solar 
corona, including, for instance, significantly higher degrees of flaring activity, and 
higher coronal plasma densities (see sec !§■ 


(a) Activity-rotation relationship 


Since observations of stars c.1980 with the X-ray observatory EINSTEIN (Gi- 


acconi et al. 1979), it became clear that rotation was the most important stellar 


parameter determining the observed activity level (Pallavicini et al. 1981). There 
is some observed residual dependence on mass of the saturation threshold, and it 
was discovered that the rotation-activity relationship becomes tighter when us¬ 
ing instead of rotation period a Rossby number Ro=P rot /T , where r is the mass- 


dependent convective turnover time (Noyes et al. 1984 Pizzolato et al. 2003). 


Wright et al. (2011) analysed the X-ray activity and rotation relation for a large 


sample including 824 late-type stars with accurately determined X-ray luminos¬ 
ity, Lx, and P ro t- The relation between Lx/L boi and Rq reveals three different 
activity regimes: unsaturated ( Rq> Ro sa . t = 0.13), saturated (i?o< Ro sat , and log( 
Lx/L boi) ~ —3.13), supersaturated (Rq< 0.01; Wright et aL| [2011). 




p rot/ da y s R ° 

Figure 2. Relationship between Lx/L boi and P ro t {left), and Lx/Lb 0 i and Ro (right)-, from 


Wright et al. (2011). For Ro > 0.1 the stellar activity scales with Ro. For Ro below 


this threshold value the activity level is essentially constant at a saturation value of log( 
Lx/L(,oi) ~ —3. 


The relation between rotation and convection parameters and non-thermal losses 
from coronal and chromospheric activity provided strong evidence that the latter 


quantities can be used as proxy indicators of the efficacy of stellar dynamos (Ros- 


ner, 

1980 

Noyes et al. 

1984 


1984). Dynamo generation of magnetic fields is commonly 


parameterized by a dynamo number (No) which is essentially the ratio between 
magnetic field generation and diffusion. The Rossby number was shown by |Noyes| 
et al. (1984) to be related to the dynamo number for an a-ui dynamo distributed 
in the convection zone according to the relation No = Rq 2 . However, Montesinos 


et al. (2001) have pointed out that in the currently favored interface dynamo acting 


at the tachocline (e.g., Parker 1993 Charbonneau & MacGregor 1997), the rela¬ 


tion between the dynamo number and the Rossby number is not quite so simple. 
They showed that the dynamo number does not depend on the Rossby number 
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alone, but also on a term including the differential rotation across the tachocline, 

AQ. 

In the unsaturated regime, the slope of the Lx/Lvs. Rq distribution then 
has potentially interesting implications for dynamo theory. The steep value of the 
slope appears to be inconsistent with the distributed dynamo but can be more 
easily accomodated by an interface dynamo acting at the tachocline. For stars of 
fairly similar spectral type, Wright et al. (20111 reduced the Montesinos et al. 
(2001) dynamo number formulation to Njj oc Afi/(fii?o). If differential rotation 
driving the dynamo is proportional to the rotation rate, Lx/L boi should vary as the 
square of the rotation rate (and similarly, with 1/Rq), as was originally observed 
by Pallavicini et al. (1981). Reiners et al. (2014) note that the approximate depen¬ 
dence of X-ray luminosity on P^ 1 / 2 for the unsaturated regime indicates that the 
convective turnover time scales approximately proportionally to L bol . However, 


Wright et al. (J2011 ) found instead a best-fit slope of 2.7. While this is somewhat 


dependent on the stellar sample, the observed slope implies a decrease of differential 
otation (AH oc fi 2 / 3 ) as the star spins down. This relation is consistent with those 
lerived from the limited differential rotation data for single dwarfs (AH ~ H°- 7 , 

it becomes problematic in 


Donahue et al. 

1996] AH ~ H 10 , 

Saar 

2011) 


the light of recent measurements of differential rotation based on high-quality Ke¬ 
pler photometry. Those data instead indicate AH to be essentially constant over a 
fairly wide range of rotation rate H (e.g. Kiiker et al. 2011 Reinhold et al. 2013). 
Either the interpretation of the Kepler photometry is incorrect, or the elementary 
dynamo theory will need to be revised to explain the lower activity end of the 
rotation-activity relation. 

We point out that the observations here are of latitudinal differential rotation, 
while the component most relevant for dynamo action is the radial component. 
The two are connected, and likely scale together, but their relationship is non¬ 
trivial (e.g., radial AH is roughly constant on cones in the Sun, but is expected 
to be constant on cylinders in faster rotators). Astroseismology gives some hope of 


resolving this issue, at least for evolved stars (Deheuvels et al. 2014); radial AH 


measurements in dwarfs may still be very challenging (e.g., Lund et al. 2014). 

The saturation value of Lx/L;, 0 ; above which faster rotation results in no further 
increase in Lx, appears independent of spectral type. The causes of the observed sat¬ 
uration of X-ray stellar activity are not well established, though several candidates 
have been identified. The main underlying question to address is whether the satu¬ 


ration of activity is caused by a saturation of the dynamo efficiency (e.g., Gilman 


1983), or is rather a separate consequence of the fast rotation, such as for instance 


saturation of the coronal filling factor (e.g., 

Vilhu 

1984 

) or centrifugal stripping 

of the corona (e.g., 

Jardine & Unruh 

1999). The lack of observed saturation in 


some chromospheric activity indicators (e.g., MgII; 

Cardini & Cassatella|;2007 

) and 

the saturation of others (e.g., Call; 

Mamajek & Hillenbrand 

2008 

Marsden et al. 


2009) at a different Ro sat value, are in contrast with the idea that the dynamo 


itself saturates. Analogously, observations of rotational modulation of X-ray emis¬ 
sion (e.g., Marino et al.. 20031 and small coronal filling factors (e.g., Testa et al. 


2004a) in saturated stars, indicates that saturation is likely not due to a saturation 


of the filling factor. The transition between saturated and unsaturated regime could 
instead be due to a change of regime of the underlying dynamo mechanism from 
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a convective (C) dynamo in fast rotators, where the stellar core and envelope are 


essentially decoupled, to an interface (I) dynamo for the slower rotators (Barnes 
200 36|g |. The similar shape, in the mass-age space, of the unsaturated-saturated 


X-ray activity regimes, and of the C —> I dynamo (see Fig. 4 of Wright et a/.||20lT I 
supports this hypothesis. As noted by Barnes (200361, late M-type stars (< 0.25 
M 0 ) are all in the C sequence. 

Recently, Blackman & Thomas (2015) explain saturation in terms of the rel¬ 
ative timescales for shear and convection: if the shear time-scale is shorter than 
the convective turnover time, convection becomes unimportant for magnetic field 
amplification. Towards later spectral types, stars have longer turnover times and 
saturate at larger values of P rot (see e.g., Fig. 6 and 7 of Wright ~ef gZ.||2011 ). This 
makes sense in terms of the toroidal shear rate being roughly proportional to ro¬ 
tation rate, but again is problematic if the rate of shear is constant for different 
rotation rates. 

Some extremely fast rotators show a decrease in X-ray activity below saturation 
levels (Randich et al. 1996| ). This supersaturation of the coronal activity is not 
observed in M-type dwarfs (e.g., Jeffries et al. 2011), but only in earlier spectral 
types. Indeed, F-type stars do not appear to go through a saturation phase, but 
instead evolve from a supersaturation phase directly into the unsaturated regime 


(Wright et al. 

2011 

). Supersaturation does not appear to occur for chromospheric 

activity 1 

Marsden et al. 

2009 

Jackson & Jeffries 

2010 

), challenging the possibility 


that supersaturation might be due to a change in the underlying dynamo. 

Possible theories put forward to explain supersaturation include: (1) “coronal 
stripping” - coronal loops become unstable due to centrifugal force, when the star 
is rotating fast enough that the Keplerian corotation radius gets very close to the 
stellar surface (e.g., Jardine 2004); (2) decrease in filling factor due to poleward 
migration of the active regions, caused by the effect of rotation on the internal 
radiative transfer in the star - “gravity darkening” effect - and strong convective 


updrafts in the outer convective envelope (e.g., St§pien et al. 20011. Wright et al. 


(2011) find that their sample is compatible with both these theories, though a 
stronger correlation appears to hold for the latter. 


(6) The End of the Main-Sequence 


Stars later than about M4 are fully convective and so the interface dynamo 
at the base of the tacliocline should not be applicable. There has been substan¬ 
tial effort aimed at understanding magnetic activity and coronal heating around 
this transition to fully-convective stars in order to probe changes in the underlying 
dynamo. There is no obvious transition in coronal properties or heating efficiency 


based on Lx/L b 0 i as the fully convective limit is traversed (Fleming et al. 1995 


Giampapa et al. 1996). Drake et al. (1996) showed that quiescent coronal emission 


from the M7 dwarf VB8 had remained at a similar level between 1979 and 1994, 
and argued for a turbulent dynamo as being the dominant field generation mecha¬ 


nism (Durney et al. 1993 Weiss 1993). Magnetic surface imaging shows the field 


geometry changes from mainly toroidal and non-axisymmetric (M0-M3) to predom¬ 
inantly poloidal and axisymmetric (late M; Donati et al. 2008; Mo rin et a6|[2008| , 
and it is not obvious how this geometry might be maintained by a purely turbulent 
dynamo. 
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A transition in coronal heating efficiency does exist, however. A more sensitive 
Chandra survey of very low mass stars and brown dwarfs by Berger et al. (2010) 
found that it occurs at masses significantly below the fully-convective limit, at a 
spectral type of about M8-9. The decrease of coronal emission for these late spectral 
types is possibly due to increased neutrality in the photosphere, which can hamper 
transport of magnetic energy through the stellar atmosphere (Mohanty et al.. 2002). 


(c) Flaring activity 


As discussed above, “nanoflares” are widely considered a viable candidate heat¬ 
ing mechanism for stellar coronae. The difficulty in detecting these small flares (see 
e.g., Testa et al. 2011 Testa & Reale 2012 Testa et al.. 2013 2014) renders it a 


challenge to establish their importance for the heating of the solar corona. Inves¬ 
tigations of solar flares and microflares indicate a power-law distribution of flaring 
events as a function of energy ( dN/dE oc E~ a - 


Lin et al. 

e-g; 

1984 

)• 

Hudson ( 


showed that the index, a, of the power-law distribution has fundamental implication 
on whether flares might significantly contribute to coronal heating: a > 2 implies 
that smaller events dominate the total energy dissipation rate and that coronal 
heating could be ascribed to smaller flares. Solar flare studies have yielded values of 
a between ~ 1.5 and ~ 2.7 (e.g., Crosby et al. 1993 Parnell & Jupp 2000 Benz & 


Krucker 20021, leaving the issue unsolved. A similar range of values for a has been 


determined for stellar observations (e.g., Audard et al. 2000 Kashyap et al. 2002). 


A recent study, refining the analysis of stellar flares, fitting the photon arrival time 
distribution with a model of stochastically produced flare emission, for a sample of 
~ 50 stellar spectra obtained with Chandra , shows that for active stars a is often 
> 2, possibly implying flare heating (Kashyap et al. in preparation). 

The hypothesis of flare heating being significant, if not dominant, in active stars 
is also supported by studies of stellar radio and X-ray emission. The radio emission 
from active stars (including RS CVn binaries and FK Com-type c oronae) s h ows a 
good correlation with their soft X-ray emission (Lxoc L^ 73±0 ' 03 ; Benz & Guedel 
e.g., 1994) over 10 orders of magnitude, from solar microflares (though not solar 


quiescent emission) to the high activity end of stellar coronae. The similar behavior 
of solar flares and active stellar coronae supports the importance of flare heating 
for “quiescent” emission of active stellar coronae (e.g., Forbrich et al. 2011). The 


radio/X-ray correlation is however found to break down for spectral types later 
than ~ M7, due to a substantial decrease of X-ray activity (Berger et al. 2010). 


The stellar observations described above suggest a scenario in which flare heating 
might be significant only for active stars. A scenario with different regimes possibly 
dominated by different heating processes is also suggested by other stellar data. For 
instance, studies of coronal densities from high-resolution X-ray spectroscopy show 
that active stellar coronae are characterized by plasma densities significantly higher 
than in the solar corona and other low activity stars (~ 10 10 cm -3 at log(T[MK]) ~ 
6.5 and ~ 10 12 cm -3 at log(T[MK]) ~ 6.8), implying remarkably compact hot 


coronal structures in active stars (e.g., Testa et al. 2004a similarly compact coronae 


are also inferred by different diagnostics, such as X-ray optical depth, and modeling 
of flares and Fe fluorescent emission, Testa et al.|2004fc| 2007&|a 2008). The derived 
filling factors increase with Fx, at both log(T[MK]) ~ 6.5 and ~ 6.8. At a surface 
X-ray flux level corresponding to the Fx typical of solar active regions, the filling 
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factor of ~ 3MK plasma approaches unity, while the filling factor of the hotter 
plasma (> 6MK) shows a change of slope, increasing more steeply with increasing 
Fx- The change of slope could suggest a more rapid increase of flaring activity once 
the stellar surface becomes covered with active regions (Filing factor = 0.1 — 1), 
possibly due to more interactions between active regions (e.g.,|Testa et al., 2004a). 


( d ) Evolution of activity with stellar age 

As discussed above (see beginning of sec. [2]) both solar and stellar studies have 
found that X-ray emission correlates well with magnetic flux. Several stellar studies 
have investigated the correlation of the average unsigned magnetic field ((|£>i|)), 
from Zeeman broadening (ZB) measurements, with rotation, and found that (|-Bi|) 
increases with stellar rotation before reaching a saturation level (e.g., Saar| |l996| 


2001 Reiners et al. \ 2009), similar to other activity indicators (see e.g., sec 2|a I. 


We note, however, that since (|Si|) = fB\ (where / is the intensity-weighted sur¬ 
face filling factor), the saturation of (|Si|) could occur in either / or B\, or both. 
Analysis by Cranmer & Saar (2011), however, suggests that the primary rotational 


dependence comes through /. 


Vidotto et al. (2014) have recently investigated in detail the relationships be¬ 


tween large-scale surface magnetic fields and stellar age, rotation and X-ray emis¬ 
sion level for 73 stars (0.1 <M*< 2.0M 0 ; 1 Myr < age < 10 Gyr). From Zee- 
man Doppler imaging (ZDI) observations they calculate the unsigned surface mag¬ 
netic field strength ((|l?v|); integrated over the surface of the star) and flux <f y 
(= AnRl(\By\)), based on the radial component of the observed surface field (since 
they are mainly interested in the field associated with the stellar wind), and com¬ 
pare their findings with previous works that used ZB (e.g., Saar, 1996} 200l). The 


dependence ((|-B v |)oc r 0 ' 655±0 045 ; 

Vidotto et al. 

20141 as the Skumanich law for 

stellar rotation vs. age (Q oc t -0 ’ 5 ; 

Skumanich 197 

2). Vidotto et al. (2014) analyze 


the relationships between (|l?v|) and Ro, as well as Lx/L ^oi vs. (|l?v|) and in both 
cases they find that they are similar (compatible within the uncertainties) to the 
dependences previously derived for (|Si|) ( Saar[ [1996} 2001 Reiners et al. |2009 


see Figure [3]). These results suggest that the ZDI (large scale) and ZB (large and 
small scale) field measurements are coupled, and that the small and large scale 
fields share the same dynamo field generation process and they might contribute 
similarly to the X-ray activity. A mor e recent analysis , howev er, suggests a steeper 
relation / oeRo~ b with 2.5 < b < 3.4 (Cranmer & Saar 2011), and hence a similar 
relationship for (\Bj\). These new slopes are in better agreement with the coronal 


Lx/L boi vs. R 0 (Wright et al. 2011) and the 4 >b vs. Lx/Lf, 0 ; (Pevtsov et al. 2003) 
relations, but imply a steeper age dependence than (|I3y|). This is not necessarily 
a problem; since (|Sy|) measures the net flux, it is reasonable to think that as the 
star becomes more active and crowded with magnetic regions, \By\ measurements 
will be more affected by vector field cancellation and thus rise more slowly with 
rotation than the average unsigned magnetic field. 

The evolution of the coronal activity of late-type stars with age has been the 
subject of extensive investigations, and more recently with particular focus on un¬ 
derstanding the impact of stellar activity on the evolution of planetary atmospheres 


(e.g., Ribas et al. 2005 Sanz-Forcada et al. 2011). Studies of X-ray/EUV/UV emis- 
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Figure 3. Correlations between average large-scale field strength (|-Bv|) (from Zeeman 
Doppler imaging), stellar parameters (age, and Ro), and activity (Lx/L;, 0 ;). (a): log(|Bv|} 
vs. log age, for non-accreting late-type stars, (b): log(|Bv|) vs. logRo; the dotted line is 
indicative of the correlation, found by |Saar (2001), between (|-Bi|) (from Zeeman broad¬ 
ening measurements) and Ro- (c): log(Lx/L(, 0 ;) vs. log(|Bv|); the dotted line shows the 
slope of the log(Lx/Lft 0 ;) vs. (|Bi|) (Saar] 2001 Wright et al.\ |2011|). (d): correlation of 
magnetic flux (log<l>v) with logRo- All figures are from Vidotto et al. |{2014). 


sion for a sample of solar analogs at different ages (~ 0.1 — 6.7 Gyr) revealed that 
the high-energy stellar emission declines following a power-law relationship with 
age, and the decline is faster at higher energies (i.e., the power-law slope decreases 


monotonically from the X-ray to the UV band (Ribas et al. 2005). The decline 


of coronal activity with age is accompanied by lower coronal temperatures and de¬ 
creased flaring rates ( Telleschi et al.\ 2005). Stelzer et al. (2013) have recenty studied 
the chromospheric, transition region and coronal emission (by analysing Ha, UV 
and X-ray data) of 159 M dwarfs (M0-M3) within 10 pc from the Sun, sampling a 
wide range of stellar ages, from ~ 0.01 Gyr to ~ 3 Gyr. Their findings show that 
the evolution of chromospheric and coronal activity of M dwarfs is similar to higher 
mass solar analogs, with the X-ray activity levels declining faster with age than 
the chromospheric UV emission. The faster decline with age of X-ray activity com¬ 
pared with the transition region and chromospheric activity can be interpreted in 
the light of the steeper dependence of the X-ray emission on magnetic flux density, 
as discussed above at the beginning of sec. [2] (e.g., Schrijver & Zwaan, .2008). 
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3. Activity cycles and grand minima 


The activity of the Sun is known to undergo cyclic variations with an 11-yr pe¬ 
riod (22-yr period for magnetic activity), as manifest in several activity indicators, 
including the number and surface area of sunspots, total solar irradiance, chromo¬ 


spheric and coronal activity, and magnetic morphology (e.g., Hudson 1988 Hath¬ 


away 2010). The first evidence of magnetic activity cycles in solar-like stars has 


been provided by the extensive chromospheric surveys carried out through Mount 


Wilson Can HK observations (Wilson 1978), which showed that activity cycles are 


a common occurrence in late-type stars (~ 60% of the Mount Wilson main-sequence 
F- to M-type stars have cyclic chromospheric activity with periods up to ~ 26 yrs; 


Saar & Brandenburg 

1999 

Olah et al. 

2009 


but otherwise clear correlations 
between activity cycles and stellar parameters have not been found. 

While cyclic chromospheric activity has been observed for a large number of 
stars, the detection of X-ray stellar activity cycles is significantly more challenging 
due to the intrinsic high variability of X-ray emission exhibited by late-type stars 
and the difficulty of continuous X-ray monitoring. To date, convincing X-ray cycles 
have been seen in only a handful of stars, such as for example 61 Cyg A (K5V; 
Hempelmann et al. 20061, HD 81809 (G5V; Favata et al. 2008), and a Cen B 
(K1V; Ayres 2009; DeWarf et al. 2010| ). These are all stars of relatively low magnetic 
activity levels. At high activity, the evidence to date suggests activity cycles are 
not present, at least in X-rays. The eclipsing RS CVn-type active binary AR Lac 
has been observed by all X-ray satellites since the Einstein observatory in 1980, 
and data obtained over the past 33 years indicate a remarkable constancy of X- 
ray luminosity with variations outside of obvious flares typically being about 15% 


and with a maximum of 45% (Drake et al. 2014). In this context, it might be 


expected that cyclic activity sets in at some activity level (age), as a star spins down. 
Recently, a short coronal and chromospheric cycle (~ 1.6 yr) has been detected in 
the ~ 600 Myr old solar-like star, t Hor (F8V; Sanz-Forcada et al. 2013). This is 
the youngest and most active star for which a cycle in chromospheric and coronal 
emission has been detected so far, and is therefore potentially important to constrain 
the dynamo properties in young solar-like stars and the onset of activity cycles. The 
boundaries between saturated and unsaturated X-ray emission delineated by |Wright| 
et al. (2011) indicate l Hor became unsaturated about 400 Myr ago; searching for 


cycles down to ages of 200 Myrs or so would be of great interest for probing the 
transition from saturated magnetic activity. 

An interesting aspect of the long term variability of stellar activity is the phe¬ 
nomenon of magnetic grand minima (MGM). For the Sun, long term records of 
sunspot numbers revealed the occurrence of occasional grand minima where some 
aspects of solar activity (such as sunspot number) have temporarily deviated from 
cyclic behavior, though other indicators (e.g., 10 Be abundance in ice cores) show 


a persistence of the magnetic cycle through MGM (Beer et al. 1998). The rela¬ 


tively short baseline of temporal monitoring of magnetic activity in stars hampers 
detection of MGM in stars. MGM have previously been identified as low-activity 
stars with chromospheric activity indicator R' HK = (Fhk — F p hot)/ Ebol ( Fhk is the 
raw calibrated HK core flux, K p hot is the photospheric component of the HK core 
flux, and F^oi is the bolometric flux; Noyes et al .| 1984) below a certain threshold 
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Figure 4. Left-. Distribution of Maunder minima candidates as a function of T e // (data from 
|Saar fe Testa|[2012| ). Right: X-ray surface flux Fx for a few candidates stars in magnetic 
grand minima and comparison stars (see Saar & Testa 2012 for details). Squares are MGM 
candidates, diamonds are stars in temporary minima, and the triangle is the non-MGM 
companion to the MGM candidate HD 179958. For stars we recently observed with the 
Chandra X-ray observatory ( Saar fc Testa| 20121 vertical dots connect fluxes calculated 
for a single assumed coronal temperature (Tx = 0.1 keV) or the best fit Tx (heavy solid). 


(log R' hk < 5.1; Henry et al. 1996). However, Wright (20041 noted the influence 


of gravity and metallicity on R' HK , complicating the definition of activity minima 
based on an R' HK threshold. 

We recently refined the MGM candidate definition criteria, in particular to re¬ 
move metallicity effects, and explored their X-ray emission properties in comparison 


with low-activity stars (Saar & Testa 2012). Taking into account the metallicity, 


a lower boundary for log R' HK is found as a function of metallicity (R'hk~ ~ 0-213 
[M/H] —5.125) for dwarfs ( Saar fe Testa| 2012). We therefore defined MGM candi¬ 


dates as dwarf stars (as determined from spectroscopically derived T e ff and gravity) 
lower boundary (log R' HK — log R ' HK min < 0.054), and with low 


close to this R' HK 



instance, higher activity stars, if they experience MGM events, could have MGM 
R'hk values larger than stars at the lowest activity levels. Assuming these criteria 
and using the sample of Wright et al. (2004), we find that about 7% of dwarfs are 


MGM candidates, with no evident dependence on metallicity. The distribution of 
the MGM candidates as a function of T e ff peaks near the solar T e f / (see left panel 
of Fig. [4]). 

We recently carried out a campaign for X-ray monitoring of a few of these 
MGM candidates and the derived X-ray fluxes are compared to other low-activity 
stars in Fig. [3] The MGM candidates have X-ray emission similar to the Sun at 
normal magnetic cycle minima, though their coronae seem cooler than the solar 
corona at minimum, perhaps suggesting a larger contribution from coronal holes 
rather than quiet sun-like regions in these extremely low-activity coronae. The X- 
ray fluxes of these low-activity stars can be compared with typical Fx of the solar 
quietest coronal features, such as quiet sun and coronal holes. For the quiet Sun 
Fx > 2 x 10 4 erg cm~ 2 s _1 (Chamberlin et al. 2009 Linsky et al. 
for coronal holes Fx is approximately in the range 5 — 10 x 10 ;i erg cm 


2014|), while 
e-g-, 


s ’ 1 
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Schmitt|[20"l2] estimated 
^-2 o-i 


Fx 


10 4 erg cm 2 


s 1 . while we estimated Fx ~ 7 x 10 3 


erg cm 2 s 1 by using temperature and emission measure estimates from Kano 


et al. 


et al. 


(2008) and synthesizing the emission using CHIANTI; Dere et al. 1997 Landi 


2012). Some of these low-activity stars can reach Fx values possibly even lower 


than the Fx of solar coronal holes cited above, but these appear to be unusually 
cool and weak (normal) cycle minima, rather than extended “grand minima” states 
(Saar & Testa| 2012). We note, however, that direct comparisons between solar 
and stellar values are often not straightforward, for several reasons including (a) 
different observing passbands of solar and stellar instruments, and (b) significant 
uncertainties in the Fx estimates above, reflecting uncertainties in the derivations 
of temperature and emission measure of coronal holes, e.g., due to scattering and 
contribution of overlying QS plasma, especially at the CH boundary (as well as 
other assumptions, e.g., on coronal abundances). 

Solar and stellar observations of activity cycles and cycle grand minima are 
crucial to constrain dynamo models, and investigate the activity conditions and 
other properties of MGM, their underlying causes and transient nature (i.e., what 
dynamo properties determine the onset and duration of MGM). Since the Sun is 
currently not in such a state, stars offer the only avenue of observational attack. 


4. Abundance anomalies 


Early spectroscopic observations of the solar corona revealed a chemical composition 
different from the underlying photosphere, in particular with a dependence on the 
element’s first ionization potential (FIP): low FIP (e.g., Mg, Fe, Si) elements appear 
2-4 times more abundant in the corona than in the photosphere, compared to high 
FIP elements (e.g., C, N, O; Meyer e.g., 1985 Feldman e.g., 1992). X-ray and 


EUV high-resolution spectroscopy showed that other stars also exhibit fractionation 
effects in their coronal composition compared with photospheric abundances (see 
e.g., review by Testap OlQq and references therein), and a FIP effect similar to the 
one observed on the Sun was observed for some low activity stars, such as e.g., a Cen 


(e.g., Drake et al. 1997). The chemical composition of the coronal plasma appears 
to depend on the stellar activity level, with active stars characterized by an inverse 
FIP effect, i.e., with coronal enhancement of high FIP elements (e.g. 


et al. 2001 Garcia-Alvarez et al. 2009 Testa 2010aI. The coronal abundance 


Brinkman 


patterns and their possible dependence on stellar parameters and activity levels are 
interesting in the context of coronal heating. The fractionation is thought to occur 
during transport to there corona from lower altitudes, as species transition from 
being neutral to ionized. As such, abundance anomalies offer the potential to probe 
the physical processes behind this transport, which are also likely related to the 
mechanism(s) of stellar atmospheric heating. 

Recently, additional studies have suggested that, whereas for high activity stars 
the coronal abundances depend on activity level, for low to moderate activity (main- 
sequence) stars (log L x < 29.1) they might have a stronger dependence on spectral 
type ( Wood fc Linsky| 2010| Wood et al. 2012 Wood fe Laming||2013 ). A change 
in abundance patterns between high activity and low/intermediate activity levels 
might reflect a transition between different heating regimes, as discussed in sec 
above, in the context of flaring activity. 
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Figure 5. For the source sample of |Wood et al. 1 2012| ) and Wood & Laming (20131, we 
plot FIP bias vs. Lx/L boi (left), vs. l/Ro (center), and vs. r c /Ro (right). 

To further explore the effect of activity level on coronal abundances we use the 
measurements by Wood et al. ( 2012| ) and Wood fc Laming] (]2013| and plot their 


measured FIP bias as a function of the activity index Lx/Lf, 0 ; and as a function 
of \/R q (left and center panel of Figure [5] respectively). These plots suggests that 
a dependence of chemical fractionation on stellar activity is present also for low to 
moderate activity stars, likely together with a spectral type dependence. This is 


Robrade et al. 2008 

) and in 

Drake 

2011 

Landi & Testa 


Landi & Testa 2015). High activity coronae do not exhibit such 
dependence but instead appear to show a “saturated” coronal Ne/O (|Drake fc 


Testa| 2005), as also discussed by Wood et al. ( 2012| ). Figure [5] (center panel) shows 
a general increase of FIP bias with logf?o up to saturation at 1/Rq ~ 10, with some 
outliers. Investigation of the outliers suggests an additional mass, and specifically, 
convection zone-related term might further improve the relationship. The right 
panel of Figure [5] explores an additional dependence on r (i.e., t/Rq)', the result 
shows reduced scatter, with only one star (Procyon) significantly discrepant. We 
note that the FIP saturation point (l/f?o ~ 10) approximately corresponds to the 
saturation threshold of coronal emission (Fig. |2j) , suggesting a possible connection 
between the phenomena: the maximum FIP bias might be limited by the maximum 
coronal activity possible. 

Efforts have been made in searching for a unified model to explain chemical frac¬ 
tionation in stellar coronae, in particular investigating the effect of ponderomotive 
forces associated with the propagation of Alfven waves through (or reflection from) 


Laming 

2004 

2009 

2012 


2012 Wood & Laming 2013). While 


these studies provide a promising start, further progress is needed, both on the 
observational side (e.g., a more thorough disentangling of dependence on different 
stellar parameters, improved determination of abundance anomalies in different 
solar coronal features), and on the theoretical modeling, in order to exploit the 
potential of coronal abundances as diagnostic of coronal heating processes. 


5. Discussion &; Open Questions 

We have attempted a very brief summary of some recent developments in the study 
of stellar coronal activity within a broad context of coronal heating. While the 


Article submitted to Royal Society 






































































Stellar Activity and Coronal Heating 


15 


general picture of stellar coronae has been greatly clarified over the last decade 
or so, as stellar studies become more detailed and instrumentation more powerful, 
almost as many new questions are raised by the new data as are answered. Here, 
we outline what we consider to be some of the more important outstanding issues, 
and the observational and theoretical avenues we need to take to begin to address 
them. 


1. What are the respective roles of interface dynamos and purely convectively- 
driven dynamos in stars ? The realization that magnetic dynamos ultimately 
give rise to coronal heating raised the question of how the dynamo behaves 
in stars without the means for rotational shear at the base of the convection 
zone—in fully convective T Tauri stars or very late M dwarfs for example. 
That there is no precipitate decline or change in coronal properties between 
fully convective stars and stars with convectively stable cores is perhaps hint¬ 
ing that the tachocline is not such a vital part of the process. Dynamo theory 
for non-solar cases would bear revisiting, especially in the light of differential 
rotation measurements (below). 


2. How do stellar activity and magnetic dynamos map to the growing number 
of observations of stellar differential rotation? The simplicity of the relation 
between unsaturated X-ray emission and the square of the Rossby number 
is beguiling, but elementary dynamo models depend on radial differential ro¬ 
tation scaling with rotation rate. Yet both theory and some observations (of 
admittedly, as noted in sec 2]a above, the latitudinal shear component) are 
converging on shear with only a very weak rotation rate dependence. Rota¬ 
tional shear must play a role in magnetic field generation, but in dynamo 
models a constant shear implies either a stronger rotational dependence of 
the convective term (a-effect), which seems unlikely, or else rotation plays 
a role not yet considered. The logical step would be, ironically, backwards, 
toward the idea in which the dynamo is more distributed through the convec¬ 
tion zone than seated at the tachocline. One problem is that the large Kepler 
sample of differential rotation measurements has not been studied to any ex¬ 
tent in X-rays, or with a technique that provides an unambiguous diagnostic 
of magnetic flux generation. They are also less well characterized in terms of 
their binarity or evolutionary status: Saar (2009 2011) argue that binary and 
evolved stsrs should be considered separately, as they appear to have differ¬ 
ent, reduced AO dependences. X-ray luminosities combined with differential 
rotation measurements for well-characterized stars will provide the necessary 
data to move forward. The Kepler K2 surveys will include well-studied stars 
in a variety of clusters (e.g., M35, M44, the Pleiades, Hyades and Praesepe) 
that have detailed X-ray data, and could provide an interesting step forward. 


3. How does flaring behavior depend on spectral type and stellar luminosity class? 
We tend to use the term “flare” somewhat loosely, and while the general 
properties of flares appear similar, their detailed characteristics are probably 
dependent on the magnetic field morphology- loop sizes, active region struc¬ 
ture and field strength—that in turn probably depends on underlying stellar 
properties such as convective vigor, differential rotation and surface gravity. 
Stars like the X-ray bright evolved Capella binary (the brightest stellar X- 
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ray source in the sky), and the nearby F6 sub-giant Procyon are notorious for 


showing little or no obvious X-ray flaring activity (e.g. Testa et al. 2012). This 


is telling us something about how their surface magnetic fields are configured 
and excited (or not). The energy dissipated by large solar flares is related to 
the magnetic potential of parent active regions (e.g. Schrijver et al. 2012). A 
systematic study of flare energies and their distributions on stars over a range 
of spectral types could provide us with a picture of surface magnetic potential 
energy and how it is shaped by stellar characteristics. Another aspect to in¬ 
vestigate is whether the energy partition in stellar flares differs from the solar 
case. The argument above provides some rationale for flares to behave differ¬ 
ently on different types of stars. Chandra and XMM-Newton have succeeded 
in detecting the Neupert effect, indicating the relationship between evapo¬ 
rated and/or heated flare plasma and accelerated particles, but as in the case 
of the Sun, detection of blueshifts expected from evaporating plasma in loops 
have proved elusive. Next generation missions (e.g. ASTRO-H, Athena ) car¬ 
rying calorimeter-based imaging spectrometers with high spectral resolution 
at high energies to observe H-like and He-like complexes of abundant metals 
such as Ca and Fe will be needed for further observational progress. At higher 
energies, only a small handful of giant stellar flares have been studied in hard 
X-rays, where non-thermal signatures are easier to disentangle from thermal 
emission. Current missions with hard X-ray capability are still too insensitive 
to study stellar flares routinely. Sensitivity at least an order of magnitude 
greater than that of NuSTAR would be needed to make progress. Since the 
majority of the energy in solar flares is released in white light or CMEs (see 
below), simultaneous white light and X-ray observations are needed determine 
stellar flare total energies. Kepler has also provided a rich vein of white light 
flare data, but again on stars with no X-ray monitoring. Simultaneous opti¬ 
cal and X-ray observations of flares are few and far between. In this context, 
X-ray campaigns coordinated with the K2 survey would be very valuable. 


4. How is the magnetic energy dissipated in stellar outer atmospheres partitioned 
between coronal heating, winds and coronal mass ejections? For the Sun, en¬ 
ergy dissipated in the outer atmosphere is divided fairly equally between coro¬ 
nal heating and wind kinetic energy, with coronal mass ejections accounting 


for 5 to 10% (depending on cycle phase) of the solar wind mass loss ( Webb 
& Howard 1994| Priest | 2014). Solar-like winds are very difficult to study on 


even the most nearby stars, and the scant studies to date suggest that stel¬ 
lar wind mass loss does scale with X-ray luminosity up to a certain activity 
level. The Wood et al. (2014) study of the 500 Myr old solar analog i r 1 UMa 


instead reveals essentially a solar mass loss rate at very high activity lev¬ 
els. This suggests a fundamentally different outer atmospheric magnetic field 
structure. The role of CME’s in mass loss on other stars is also very unclear. 
Scaling solar flare-CME relations to the flaring rates of stars suggests pro¬ 
hibitively large CME mass loss rates requiring unrealistically high amounts of 


energy (Drake et al. 2013). The implication is that large stellar flares put a 


much lower fraction of their energy into CMEs. A much larger survey of stel¬ 
lar winds is needed to better understand their dependence on stellar activity 
and spectral type. Unfortunately, the technique, based on Lya observations, 
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is extremely challenging and it is probably infeasible to expand the existing 
Hubble Space Telescope sample beyond a very small number of additional tar¬ 
gets. CMEs on stars cannot presently be observed directly. Their cumulative 
effect should be seen in the same observations used to probe stellar winds, 
however. The indication of much lower mass loss rates in the most active stars 
suggest a very interesting change in the way flares work, with the implication 


that the flare-CME relation is indeed fundamentally different (Drake et al. 


2013). More observations of stellar mass loss are needed for active stars to 


make progress; unfortunately this would seem to require new observational 
techniques, and perhaps new missions. 


5. How does coronal heating work near the end of the main-sequence, where 
stellar photospheres are essentially completely neutral? The observation of a 
break in the ratio of X-ray to bolometric energy output of stars at the very 
end of the main sequence is very interesting for coronal heating and provides a 
parameter space that, if we could study it in detail, might teach us much about 
how different candidate processes work. The data suggest that the amount of 
magnetic energy that goes into particle acceleration—electrons at least, that 
give rise to the observed radio signatures—does not change appreciably going 
to these very cool dwarfs, but that plasma heating is efficiently quenched. It is 
expected that heating should decline as the atmosphere becomes more neutral 
due to large resistivities and efficient field diffusion, and consequently less 
effective driving of magnetic stresses. But why should particle acceleration 
remain efficient, and isolated flaring activity be present? Much longer and 
deeper observations of nearby objects are needed to better understand their 
X-ray activity, while a wider survey covering a larger number of object would 
probe coronal heating over a wider range of parameters and rotation rates. 
Such studies are partially within reach of current missions, although would 
require expensive commitments of time. 


6. How do coronal abundance anomalies correlate with underlying stellar proper¬ 
ties and with coronal properties such as magnetic scale height, plasma temper¬ 
ature and flaring rate? Theoretical basis for understanding the FIP Effect has 
historically not been lacking—it has been easy to point to suitable regions in 
the chromosphere where the separation is likely to take place and take advan¬ 
tage of those conditions to effect it—but models have tended to fall short, and 
the avalanche of new stellar data at the turn of the millennium from Chandra 
and XMM-Newton tended to fall outside of the scenarios tailored for the Sun. 
The ponderomotive force-based model proposed by Laming and collaborators 
offers an appealing solution to the problem that can in principle explain both 
enhancements and decrements in species under different conditions. While 
the observational database is presently fairly large, the solar case illustrates 
that different regions of the same star can exhibit fairly dramatic differences 
in chemical composition. Stars are likely to do the same things, and what we 
observe is an emissivity-weighted melange of that behavior. Abundance esti¬ 
mates are often fairly uncertain and carry considerable scatter, though they 
are also likely to be time-dependent according to flares and dominant active 
region evolution and rotational modulation. Future X-ray instruments with 
larger effective areas will be able to probe the time-dependence of abundance 
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anomalies, and begin to separate the ensemble behavior into its constituent 
parts. Long commitments of observing time on present day instruments could 
also tease out secular trends, but will be difficult to justify to allocation com¬ 
mittees. Theoretical advances through more detailed and realistic modeling 
can only help the arguments to secure the observational data. 

We look forward to the next few years to see how the many facets of the coronal 
heating problem evolve, and to the new array of questions and problems that both 
theoretical and observational progress are inevitably going to uncover! 
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